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The scalar f xc and tensor f xc exchange-correlation (xc) kernels are key ingredients of 
the time-dependent density functional theory and the time-dependent current den- 
sity functional theory, respectively. We derive a comparatively simple relation be- 
tween these two kernels under the assumption that the dynamic xc can be considered 
"weak" . We expect our formula to serve as a convenient bridge between the scalar f xc 
which directly enters many applications and the tensor f xc which, due to its locality 
in space, is much easier to approximate. 
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I. INTRODUCTION 



The concept of the dynamic xc kernel f xc plays a central role in time-dependent density 
functional theory in the linear-response regime. 1 In the general non- homogeneous case, it is 
defined as the kernel of the integral transformation 



which relates the change in the dynamic xc potential 8V XC to the change in the particle 
density 5n, where u is the frequency. The knowledge of f xc allows one to obtain the density- 
response function x of interacting electrons through the relation 1 



2 



X (r, r', u) = Xks( t > r '' w ) ~ 1 7\ ~ fxc(r, r', u 



r— r 



where xks is the density-response function of Kohn-Sham (KS) non-interacting electrons, 
and e is the charge of the electron. The scalar f xc is also an important quantity in appli- 
cations, e.g., it directly enters the formulas for the many-body contribution to the stopping 
power of materials for slow ions^ and the formula for the impurity resistivity of metals.— 

It is known, however, that at finite frequency f xc (r,r' ,u) is a strongly nonlocal function 
of its spacial variables.- This non-locality hinders the construction of easy-to-use approxi- 
mations in time- dependent DFT. The need to overcome this difficulty has stimulated the 
development of the time-dependent current-density functional theory (TDCDFT).- The key 
quantity of TDCDFT is the tensor xc kernel f xc defined as 



SA i>xc (r, u) = J fik,xc(r, r', u) Sj k (r', u) dr', 



where 5A XC is the change in the xc vector potential and <5j is the change in the current- density. 
It has been established that the tensor xc kernel f xc (r,r' ,u) is a much more local function 
of its space variables than its scalar counterpart is. As a result, the LDA to TDCDFT has 
been developed.-^ 

It would be natural to expect that any local or semi-local approximation to f xc (see, for 
instance, Ref. Q) can be translated to a corresponding approximation for the scalar f xc , 
which is more immediately useful in many applications. The relation between f xc and f xc 
is, however, nontrivial, and it has been established in the general case only recently. In its 
exact form it reads^ 

U = - ^ v- 2 v • |/, c + (xks - L) [f (xk S - L) f' 1 




Here xks is the KS current-density response function (a tensor), T is the projector on the 
sub-space of transverse vectors, V is the gradient operator, V -2 is the inverse Laplacian, and 
c is the velocity of light. It must be pointed out that the presence of the KS current-density 
response function in this expression is unavoidable as long as there is more than one spatial 
dimension (for the one-dimensional case, see Eq. (fTUl) below) .- 

Although Eq. (CQ) is exact and general, its direct use is not easy, since it involves repeated 
inversions of tensor integral operators. The purpose of this Communication is to propose 
a simpler and more practical relation even at the price of imposing some restrictions. As 
it will be shown below, Eq. ([TJ) can be, indeed, considerably simplified by the expansion 
of the right-hand side to the first order in f xc . Physically, this amounts to the intuitively 
clear assumption that f xc is "small" compared with the KS contribution, the latter entering 
Eq. (DO) through SCks- This assumption also underlies one of the most successful approaches 
to the calculation of excitation energies, namely the perturbative treatment of f xc in the 
single-pole approximation of Gross, Dobson, and Petersilka.-^ Our final result reads 

C * 

fxc = Xks^ ■ Xks • fxc • Xks • ^Xks- ( 2 ) 

We find it worthwhile to put Eq. (j2J) in words: It says that to apply f xc to a scalar 
function, first a scalar operator Xks mus t he applied, then the gradient of the result is 
taken, then a tensor operator xks is applied to the vector result of the previous operation, 
then a tensor operator f xc and xks are consecutively applied, again producing a vector, then 
a divergence of this vector is found producing a scalar, and, finally, a scalar operator x~ks 
is applied. We hope that the availability of this relatively simple expression will open the 
way to the practical use of more sophisticated xc kernels in standard applications of time 
dependent DFT. 

In Sec. HIl we provide the derivation of Eq. (j2J). In Sec. IHIl we discuss the result and 
draw conclusions. Explicit expressions for KS response functions, in the form convenient to 
use with our result, are collected in the Appendix. 
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II. DERIVATION OF EQ. fl2J 

To prove Eq. (j2J) we will need the following 

Lemma. For the response function x the relation holds 

(fx- 1 ?) 1 = x + — 2 x ■ Vx" 1 v • x- 

V / ear 

The same relation holds for the KS response function Xks- 

Proof. We can write 

X = LxL + Lxf + fxL + fxf, 
where L is the longitudinal projector 

L tj = V- 2 V,V r 

The same expansion can be performed for the inverse operator x~ X ■ Then 

i = xT 1 



= LxLx'L + LxLx'T+LxTx'L+LxTx'T 

+ fxLx- 1 L+f X L X - 1 f+fxfx~ 1 L+fxfx~ 1 f- 
Multiplying this by T from the right, we have 

f = LxLr'f+Lxfr'f+fxLr'f+fxfx- 1 ?. 

Multiplying the last equality from the left by T and by L, we will have, respectively, 

f = TxLx'T+fxfx 1 ?. 

and 

6 = LxLr'f+Lxfr'f. 
Combining Eqs. (jSJ) and (JBJ) yields 



TxT-TxL[LxL) LyT 



-i 



Tx X T = T, 



and therefore 



The operator in the square brackets in Eq. (J7J) is purely transverse, as one can verify by 
applying the L operator to its left and to its right and obtaining zero. We can, therefore, 
drop the T operators in the right-hand side of Eq. (J7J) , which leads to 

(fx-'f) 1 = X ~ X-L (LxL) ~* Lx. (8) 

We now recall the relation between the scalar and tensor response functions 

X = ^V-x-V, (9) 

ear 

which, combined with Eqs. (IE]) and (@|, yields Eq. ([3]). □ 
Now we expand the right-hand side of Eq. ([I]) to the first order in f xc : 

xL + Xks (TXksT)' 1 L (fxKsfy'x^yVV- 2 . 
Equation ([2]) immediately follows from the above equation by use of Eq. ([3|). 



III. DISCUSSION AND CONCLUSIONS 

We point out that similar to its full form of Eq. (CE]), in the simplified form of Eq. (TSJ), 
the scalar f xc depends not on the tensor f xc only, but also on the KS response function. 
However, in the case of purely longitudinal f xc (e.g., for ID inhomogeneity) Eq. (T5]) reduces 
to 

L = —^V/xoV. (10) 

eu 2 



This agrees with the expression suggested in Ref. [11|. We, however, note that in the inho- 
mogeneous 2D and 3D cases, there is no reason for f xc to be purely longitudinal (and it, 
indeed, explicitly is not such within LDA-). 

It is known that at isolated frequencies the operator xks can have zero eigenvalues and, 
therefore, be non-invertible.— ^ This poses an interesting question whether or not the scalar 
f xc can have a singularity due to the presence of Xks in Eq- © even if the tensor f xc is non- 
singular. This evidently does not happen in ID case, when the presence of xks in Eq- 
compensates the possible singularity leading to Eq. (jTOl) . However, in the general case, we 
cannot rule out the possibility that a scalar potential that produces no density response in 



the KS system (i.e., a null eigenvector of Xks) m &y nevertheless produce a finite transverse 
current response in the same system. If this happens, then the singularity in Xks remains 
uncompensated, and f xc may have singularities that are not present in f xc . This kind of 
singularity is, however, impossible at complex frequencies with finite imaginary part, since 
in this case the density-response function xks is invert able.— 

In conclusion, we have considerably simplified the relation between the two key quantities 
of the time-dependent density functional theory and the time-dependent current density 
functional theory: the scalar f xc and the tensor f xc , respectively. This has been achieved 
at the price of assuming the dynamic exchange-correlations to be weak on the background 
of the Kohn-Sham response, and solving the problem to the first order with respect to 
the former. We hope that the availability of this approximation, similar in spirit to the 
perturbative approximation in the Gross-Dobson-Petersilka approach, will stimulate the use 
of more accurate xc kernels in standard applications of time-dependent DFT. 
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Appendix 

Apart from the tensor xc kernel f xc , which is considered an input quantity within the 
context of this paper, to use Eq. (jSJ) one needs the tensor xks and the scalar xks- While 
the explicit forms of the latter two operators are well known, the purpose of this appendix 
is to conveniently represent the construct Xks • V, which enters Eq. (j2J). The tensor KS 
response functions can be written as 



XKS,ij(r,r',uj) = — n (r)5(r - r') c% - -x 

cm Acm z 

CO Eq ~\~ dry i vT\ 
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where vf) a (r) and e a are KS wave-function and eigenenergy, respectively, in the state a, f a is 
the occupation number of this state, and r\ is an infinitesimal positive. Applying V operator 
to Eq. (lA.ip from its right, gives us a convenient form of the operator Xks ' V which enters 
Eq. © 

XKStii (r, r , U ) V, = — £ 

x «(r)V^(r) - ^(r)V,C(r)] ^(r')^(r'). (A.2) 

The fact that ifj a (r) satisfy Schrodinger's equation with the eigenvalues e a has been used. It 
can be easily verified that applying V to Eq. (|A.2[) from its left and using Eq. (j9J) leads to 
a known expression for the KS density-response function 

^u-ep + e a + iri 
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